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We explore the influence of geometry variations on the structure and the time-dependence of the mag- 
netic field that is induced by kinematic a 2 dynamos in a finite cylinder. The dynamo action is due to an 
anisotropic a effect which can be derived from an underlying columnar flow. The investigated geometry 
variations concern, in particular, the aspect ratio of height to radius of the cylinder, and the thickness of 
the annular space to which the columnar flow is restricted. Motivated by the quest for laboratory dynamos 
which exhibit Earth-like features, we start with modifications of the Karlsruhe dynamo facility. Its dynamo 
action is reasonably described by an a 2 mechanism with anisotropic a tensor. We find a critical aspect ratio 
below which the dominant magnetic field structure changes from an equatorial dipole to an axial dipole. 
Similar results are found for a 2 dynamos working in an annular space when a radial dependence of a is 
assumed. Finally, we study the effect of varying aspect ratios of dynamos with an a tensor depending both 
on radial and axial coordinates. In this case only dominant equatorial dipoles are found and most of the 
solutions are oscillatory, contrary to all previous cases where the resulting fields are steady. 
Keywords: Dynamo; a effect; magnetic field orientation 

1 Introduction 

It is generally assumed that columnar flows in the Earth's outer core play an essential role 

for the generation of the geomagnetic field. At the surface of the Earth, the magnetic field 
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structure has almost an axial dipole (AD) structure closely aligned with the Earth's rota- 
tion axis. Direct numerical simulations of the geodynamo have successfully reproduced 
many observed features like, e.g., the dominance of the axial dipole and the occurrence of 
reversals (e.g. Olson et al. 1999, Ishihara and Kida 2002, Aubert and Wicht 2004, Wicht 
and Olson 2004 and references therein). The poloidal part of the field is thought to be 
produced from the toroidal part by the a-effect generated by the columnar flows, while 
the toroidal component of the Earth's magnetic field is associated to the fi-effect, but also 
again to an a-effect or even to both mechanisms together. These types of magnetic field 
generation are usually referred to as afl, a 2 and a 2 VL dynamos, respectively. 

It was one of the motivations of the Karlsruhe dynamo experiment to study an Earth- 
like magnetic field generation process in the laboratory (Gailitis 1967, Busse 1975, Stieglitz 
and Miiller 2001). However, in contrast to the axial dipole (AD) of the Earth, the eigen- 
field structure of the Karlsruhe dynamo is an equatorial dipole (ED) what has been pre- 
dicted in terms of the mean-field theory with an anisotropic a effect (Radler et al. 1998). 
Actually, a general tendency of anisotropic a 2 dynamos to produce fields with dominant 
equatorial dipole structure has been known for long (Radler 1975, Radler 1980, Riidiger 
1980, Riidiger and Elstner 1994). 

It is also well known that a transition from equatorial to axial dipoles can occur if 
some differential rotation is added (Radler 1986, Gubbins et al. 2000). However, an axial 
field orientation can also result from a 2 dynamos if the magnetic diffusion is enhanced by 
small scales of the flow (Tilgner 2004). 

Besides the axial and equatorial dipole, the quadrupole structure seems to play also 
a certain role in geodynamo models. In many kinematic models one finds a quasi- 
degeneration with the dipole field (Gubbins et al. 2000). This degeneration is also re- 
sponsible for the appearance of hemispherical dynamos in dynamically coupled models 
(Grote and Busse 2000). In this case, both quadrupolar and dipolar components contribute 
nearly equal magnetic energy so that their contributions cancel in one hemisphere and add 
to each other in the opposite hemisphere. The interplay between the nearly degenerated 
(Gubbins et al. 2000) axial dipole, equatorial dipole, and quadrupole was used in various 
models to explain the reversal phenomenon of the geodynamo (Melbourne et al. 2001). 
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With the same focus on field reversals, the importance of transitions between steady 
and oscillatory solutions of kinematic dynamos has been highlighted by several authors 
(Weisshaar 1982, Yoshimura et al. 1984, Sarson and Jones 1999, Phillips 1993, Riidiger 
et al. 2003). In an extremely reduced reversal model dealing only with the axial dipole 
it was shown that many features of reversals (typical time scales, asymmetry between 
slow dipole decay and fast recovery, bimodal field distribution) can be understood by the 
magnetic field dynamics in the vicinity of transition points between steady and oscillatory 
solutions (Stefani and Gerbeth 2005, Stefani et al. 2006a, Stefani et al. 2006b). The 
main ingredient of this reversal model, as well as of the reversal model of Giesecke et al. 
(Giesecke et al. 2005a), is a sign change of a along the radius which brings into play a 
coupling between the first two radial eigenfunctions of the axial dipole field. It should be 
noticed that such a sign change results indeed from simulations of magneto-convection 
(Giesecke et al. 2005b). 

With this background, we investigate in the present paper various kinematic dynamo 
models within cylindrical geometry. Our focus will lay first on the dominance of field 
structure: equatorial (ED) or axial dipoles (AD) or even quadrupoles (Q), and second on 
the occurrence of oscillatory solutions. The cylindrical geometry, which might seem awk- 
ward from the purely geodynamo perspective, is quite natural from an experimentalist's 
viewpoint. One could ask, e.g., how the geometry and the arrangement of spin-generators 
in the Karlsruhe dynamo could be modified in order to make its eigenfield prone to rever- 
sals. 

After presenting the general framework, we will explore geometrical effects that could 
lead to dominant AD fields in cylindrical anisotropic a 2 dynamos. The utilised numerical 
code, which is based on the integral equation approach to kinematic dynamos (Stefani 
et al. 2000, Xu et al. 2004a, Xu et al. 2004b, Xu et al. 2006), was already used for 
the simulation of various cylindrical dynamos, including the VKS dynamo experiment in 
Cadarache (Stefani et al. 2006c). 

The geometrical variations which are actually considered are the aspect ratio of height 
to radius of the cylinder and the width of the annular space to which the dynamo source 
is restricted. First we consider the geometry of the Karlsruhe dynamo experiment. Its 
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steady dynamo field is generated by a bundle of axially invariant helical columns, which 
is well described within mean-field theory as an a 2 dynamo with anisotropic a-effect. 
We find that for this experiment a dominant AD field could be achieved below a critical 
value of the aspect ratio which is not so far from the one of the real facility. In a next 
step, we explore more complex structures of a which have been derived from a flow 
described by axially invariant helical columns which are restricted to an annular space. 
The resulting a coefficients acquire a radial profile which depends on the flow structure. 
As for the modified Karlsruhe case, dynamo solutions show dominant steady AD fields 
below a critical value of the aspect ratio. In contrast to this, the reduction of the thickness 
of the annular space does not lead to a transition from non-axisymmetric to axisymmetric 
modes, although the critical dynamo numbers for both modes seem to converge. Finally, 
we have considered axial-radial dependence of a. The dynamo action works in a fixed 
annular space and again the aspect ratio of height to radius of cylinder is the varying 
geometrical parameter. In this case, non-axisymmetric oscillatory fields are the dominant 
solutions. 



2 The general concept 

We consider an incompressible steadily moving fluid with velocity u, which is confined to 
a cylinder and surrounded by vacuum. The fluid has homogenous electrical conductivity a 
and magnetic permeability \x. The fluid motion induces a magnetic field B which extends 
in whole space. The magnetic field is governed by the induction equation 

r]V 2 B + V x (u x B) -d t B = 0, VB = 0, (1) 

where rj is the magnetic diffusivity defined by rj = 1/ fia. In the mean field approach, each 
quantity is decomposed into a mean part (denoted by an overline) and a fluctuating part 
(denoted by a prime). Referring to a cylindrical coordinate system (s, </?, z), we define 
mean fields by averaging over (p. 

As we are only interested in the induction effects originated by the fluctuating part u' 
we assume that the mean motion u is equal to zero. In this case, the mean part of the 



5 



induction equation (1) reduces to 



??V 2 B + V x £ - <9 t B = , V • B = , 



(2) 



where £ = u' x B', is the mean electromotive force (e.m.f.) which is the source of 
generation of the large scale magnetic field B. This e.m.f. results from the interaction of 
motion and magnetic field at small scales. 

2.1 Representations of £ 

We consider different forms of £ which are generated by flows organised in columnar 
vortices parallel to the vertical axis of the cylinder. In the following only the a effect 
that results from these columnar structures is considered as the main contribution to the 
generation of £, other effects are just neglected. In a strict sense, such a reduction of £ to 
an a-effect term is only possible if the spatial variations of B are sufficiently weak. 

We consider first the mean e.m.f £ produced by the flow in the Karlsruhe dynamo 
experiment (Radler et al. 1998 ). Its most simplified analytical representation is given by 



where a is constant in the cylindrical volume and e z is the unit vector in axial direction. 
We point out the anisotropy of the a-effect as represented in (3). 

The next considered example of £ results from an axially invariant flow organised in 
columnar vortices equally distributed in an annular region. Similar flow structures were 
recently discussed in the context of quasi-geostrophic dynamos (Schaeffer and Cardin 
2006). A detailed description of such "rings of rolls" and the a tensor resulting from them 
has been derived by Avalos et al. (2007) and is given in Appendix A. The mean e.m.f. £ 
produced by such a flow is given, under the assumptions mentioned in Appendices A and 



£ 




(3) 



B,by 



(4) 



with the subscripts k and A standing for s, ip, or z. 
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For some flow configurations, it has been shown (Avalos et al. 2007) that the resulting 
matrix a K \ is of the form 



a ss (s) 





« w (s) 



a zip (s) 



(5) 



We have also considered an additional axial dependence of the components in (5) 
multiplying them by harmonic functions of z that vanish at the top and the bottom of the 
cylinder. This was motivated by the fact that for rolls in real rotating bodies a North-South 
antisymmetry of the axial velocity is expected, while the horizontal velocity components 
are expected to be symmetric with respect to the equator. Admittedly, the correct treat- 
ment of this problem would require a new derivation of the a matrix for such rolls along 
the lines outlined in the appendices. As a sort of compromise we focus here only on the 
general symmetry properties of the elements of the a matrix. Since a ss and a w depend 
on products of axial and horizontal velocity components, we expect an antisymmetric be- 
haviour. On the other hand, a zip should remain North-South symmetric since it depends 
on horizontal velocity components only. 

The cylinder is assumed to extend over the axial interval —H/2 < z < H/2. If we as- 
sume u s and u v to be proportional to cos(nz/H) and u z to be proportional to sm(27rz/H), 
then the new a matrix is given by 



a ss (s) cos(nz/H) sm(2nz/H) 

a w (s) cos(nz/ H) sm(2,7:z/ H) 

a zip (s) cos 2 (nz/H) 



• (6) 



Evidently, the resulting diagonal elements of a KX are anti-symmetric with respect to z = 
0, whereas the non-diagonal element is symmetric with respect to z = 0. 

Though the representation of £ defined above was derived for an infinitely extended 
conducting fluid, we assume that it applies also to a finite cylinder. This approximation 
has been successfully used, e.g. in Radler et al. (1998, 2002), to solve the Karlsruhe 
dynamo problem in a finite cylinder. In this case the symmetry of the most easily excited 
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magnetic field mode was found to be independent of the conductivity outside the cylinder, 
while the other properties of this mode well depend on the conductivity in outer space. 

3 Dynamo solutions 

Once we have defined different representations of £, we solve the mean-field dynamo 
problem in a finite cylinder enclosed by vacuum using a numerical code based on the 
integral equation approach (Stefani et al. 2000, Xu et al. 2004a, Xu et al. 2004b, Xu 
et al. 2006). The magnetic field is determined by a self-consistent solution of the Biot- 
Savart equation together with a surface integral equation for the electric potential at the 
vacuum boundaries. For time-dependent solutions, the model has to be completed with 
an integral equation for the magnetic vector potential. All field quantities are expanded 
in harmonic modes (~ exp (irrup)) in azimuthal direction and vary in time t according 
to exp (pt) with a constant p that is, in general, complex. Then, there are two ways to 
solve the integral equation system. For steady eigenfields (i.e. marginal eigenfields which 
are non-oscillatory) it is treated as an eigenvalue equation for the critical value of a. For 
unsteady eigenfields (including marginal eigenfields which are oscillatory) the integral 
equation system is treated as an eigenvalue problem in p: dynamo solutions corresponding 
to exponentially growing magnetic fields are characterised by a positive real part of p . In 
Appendix C more details about this numerical approach are given. 

We stress that the a effect has been determined under the assumption of an axisym- 
metric mean magnetic field with m = 0. Using the same a effect for other m modes is an 
approximation which is valid only if m <C n, where n is the number of pairs of rolls. In 
that case the azimuthal variation of B is weak compared to the one of u. 

3.1 Karlsruhe geometry 

It is well known that the main generation mechanism of the Karlsruhe dynamo experiment 
is an a effect which maintains, in the marginal case, a steady equatorial dipole (ED) field, 
i.e. a mode with m — 1. For numerical studies, a simplified geometry has been assumed 
in form of a finite cylinder with height H and radius R. We use this simplified geometry 
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and the £ given by (3) to compute dynamo solutions for different ratios H/ R. In figure 
1 we represent the threshold C c a of the dynamo number C a = fxaRa corresponding to 
3?{p} = 0. This is done for the two leading axisymmetric modes with m — 0, i.e. for 
the axial dipole (AD) and the quadrupole (Q), as well as for the first non-axisymmetric 
mode (m = 1) which represent an equatorial dipole (ED). All these modes are steady at 
the marginal point. 
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Figure 1: Critical dynamo number C c a as a function of the aspect ratio H/R for a 
Karlsruhe-type dynamo. The two first axisymmetric (m = 0) eigenmodes AD (axial 
dipole) and Q (quadrupole) are compared with the first non-axisymmetric (m = 1) eigen- 
mode ED (equatorial dipole). 

We have found a critical aspect ratio H/R = 0.75 which distinguishes between dom- 
inant ED and AD fields. Above this critical value ED fields are dominant while below 
this value AD fields are dominant. Actually, the critical aspect ratio of 0.75 is not very far 
from the experimental one, which is 0.83. 

3.2 Ring of rolls 

In the following, we investigate anisotropic a 2 dynamos in an annular space defined by a 
gap width of 25 R with 5 < 1. Note that in the following R refers to the radius in the mid- 
dle of the gap, and not to the outer radius. We consider both the ^-independent case with 
a given by (5) and the ^-dependent case with a given by (6). In each case we have consid- 
ered two types of flow distinguished by the radial dependence of their vertical velocities 
as defined in Appendix A. 2. We have called them FW1 and FW2. We computed the crit- 
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ical value C c a of the dynamo number C a = \10Ra where a stands for J (cn^J with (■ ■ ■) 
understood as an average over s. According to Avalos et al. (2007) the relation between 
a and the real velocity of the flow is given, under the first order smoothing approxima- 
tion, by a ps (r]5 / 'R)R Tn ±R m \\. The quantities R m ± = u ±R/r] and R m \\ = u \\R/r] are 
the magnetic Reynolds numbers expressed in terms of the characteristic velocities in the 
horizontal (i.e. perpendicular to the z-axis) and in the axial (i.e. parallel to the ,2-axis) 
direction, respectively. 

3.2.1 z- independent case 

In figure 2 the dynamo threshold is plotted in dependence on H/R for both flows FW1 
and FW2 for § = 0.5 and n = 4. Quite similar to the Karlsruhe dynamo case, a critical 
aspect ratio Hj R is also found here for both flows, which distinguishes between domi- 
nant ED and AD fields. Another critical value of H/R is found where the second (i.e. 
subdominant) eigenmode is switching between AD and Q. 

For a given ratio H/ R we found that the dynamo threshold increases monotonically 
when we reduced the magnitudes of a ss and a w while keeping a zip unchanged. This is 
related to the impossibility of having dynamo action with a z-independent horizontal flow 
only. 




H/R H/R 

(a) (b) 

Figure 2: Critical dynamo number C c a as a function of the aspect ratio H/R for an a 
matrix according to (5) with 5 = 0.5. The two first axisymmetric (m = 0) eigenmodes AD 
(axial dipole) and Q (quadrupole) are compared with the first non-axisymmetric (m = 1) 
eigenmode ED (equatorial dipole). Plot (a) for FW1 and (b) for FW2. 
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In figure 3, the rescaled dynamo threshold bC c a is plotted in dependence on 5 for both 
flows FW1 and FW2. We introduce here another distinction between the case of "free 
rolls" (for which the number of pairs of rolls is kept equal to 4 independently of the value 
of S) and the case of "compact rolls" (for which the rolls have the same extension in 
azimuthal and radial direction and the number of pairs of rolls scales like n = ir/28). In 
neither case was there any indication for a critical value of 5 below which the dominant 
m — 1 mode is clearly replaced by a dominant m = mode. However, for small values 
of 8, the values of 8C c a for the m = mode come very close to those of the m — 1 mode. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

8 8 
(a) (b) 



Figure 3: The rescaled dynamo threshold 8 in dependence on 8 for FW1 (upper curves) 
and FW2(lower curves), and fix aspect ratio H/R = 2. The dashed (solid) line corre- 
sponds to axisymmetric (non axisymmetric) fields. Plot (a) for the case of "free rolls" and 
(b) for the case of "compact rolls". 

In the case of free rolls it is remarkable that 8C c a decreases with 8 for FW1 and in- 
creases for FW2. The geometries of the magnetic field produced by FW1 and FW2 for 
8 = 0.3 have indeed different symmetries. This is illustrated in figures 4 and 5 where 
poloidal vectors and azimuthal contour of the magnetic field are plotted. On the other 
hand the symmetries are similar for 8 = 0.9. 

3.2.2 ^-dependent case 

In figure 6, C c a is plotted in dependence on H/R in the z dependent case for 8 = 0.5 
and n = 4. We find now that oscillatory non- axisymmetric fields are always the most 
easily excitable solutions for both flow types FW1 and FW2. However, the axisymmetric 
solutions are getting closer to non- axisymmetric ones as H/R is reduced. For the FW1 
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s/R 



s/R 



0.3 




Figure 4: Poloidal field component together with contour plots of azimuthal field for FW1 
with free rolls, for S = 0.3, 5 = 0.9, m — 0, m — 1. Ordinate axis runs from —H/2R 
to H/2R, corresponding to H/R = 2. For the non-axisymmetric case (m = 1) the 
plot represents the numerical solution in the meridional plane at (p = 0. However, one 
should notice the degeneration of the eigenvalue problem with respect to any rotation of 
the eigenmode in (^-direction. 

flow a transition between steady and oscillatory magnetic fields is observed for the mode 
m = at a certain value H/R ~ 1 (the precise value could not be determined since 
the numerical solution of the problem is quite time consuming) . In all cases non-dipolar 
fields only were found. 

4 Conclusions 

We have explored the influence of geometrical parameters on spatial structure and tempo- 
ral variations of magnetic fields generated by kinematic anisotropic a 2 dynamos working 
in a finite cylinder. The a coefficients were calculated for specific flow patterns, following 
the lines of mean field concept, and the corresponding dynamo solutions were calculated 
using the integral equation approach. The obtained results show that this kind of dynamos 
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s/R 



s/R 




Figure 5: Same caption as figure 4 but for FW2. 

can switch from dominant equatorial dipoles to dominant axial dipoles just by reducing 
the aspect ratio of the cylinder. This transition occurs for quite different forms of a: con- 
stant, as in the Karlsruhe dynamo experiment, or having a purely radial dependence, as 
the one obtained in a flow described by axially invariant helical columns. On the other 
hand, such a transition does not occur when the relative gap width 5 is reduced (at least 
not for the considered aspect ratio). When a has an additional axial dependence, domi- 
nant dynamo solutions are only oscillatory m—1 modes. In addition for the m = mode 
both steady and oscillatory solutions were obtained. 
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Il/K H/R 

(a) (b) 

Figure 6: Critical dynamo number C„ as a function of cylinder aspect ratio H/R for an 
a matrix according to (6). (a) FW1 and (b) FW2. The dashed (solid) line corresponds 
to axisymmetric (non axisymmetric) fields. Solutions in (a) are always oscillatory except 
for m = and higher values of H = 1.04 where solutions are steady. Solutions in (b) are 
always oscillatory. 

Appendix A: Specification of the velocity field 
A.l General assumptions 

We specify the motion of an incompressible conducting fluid u, so that it corresponds to 
a ring of columnar vortices. The ring is coated by an interval 1 — 5<s/R<l + 5 with 
5 < 1. Outside this interval the fluid is assumed to be at rest. It is assumed that u is 
steady, z-independent and varies with ip like exp(in^), where n is the number of vortex 
pairs. We use the representation 

u = -V x (e 2 x V$) -e 2 x V\&, 

$ = m || R 2 (fi(s) cos(ntp) , ^ = Uq±R if)(s) cos(ny). (A.l) 

The two terms on the right-hand side of u correspond to the vertical (poloidal) and hor- 
izontal (toroidal) parts of the velocity. The constant quantities uq± and uq\\ define the 
intensity of the considered flow. We further express u by 

u s = u s (s) sin(nyj), = u v (s) cos(rup), u z = u z (s) cos(rup). (A. 2) 
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The connection between (A.l) and (A. 2) is given by 



Tb dip 
u s = -u ±R -ip , Uip = -u _i_R t^- , u z = -u Q \\R 2 D n (f) , (A.3) 



where D n (f) = s 1 d s (s d s 4>) — (n/s) 2 4>. 

A.2 Specific examples 

We consider two flows which differ only in the radial dependence of u z . The first flow 
(FW1) is defined by 

^ = c^(i-e)\ uju H = c z (i-e)\ z= (s/R) 6 ~ l , if iei<i 

= ^ = otherwise . (A. 4) 

The second one (FW2) by 

^ = c,{i-e)\ ^ = c,(i-e)\ n = { -^l zl , if ki<i 

= ^ = otherwise . (A. 5) 

The factors , and C z were chosen such that the average of u z /u \\ over a surface 
given by 1 — 5 < s/R < 1 + 5, —ii/2n < (p < ir/2n and z/R = constant as well as the 
average of u^/uq^ at y? = over 1 < s/R < 1 + 5 are equal to unity, 



15n5 7 



25 (15 - 405 2 + 33<5 4 ) + 15 (l - 5 2 f log (^-^ 



(A.6) 



The flow definitions (A. 4) and (A. 5) ensure that u s , u v and u z are continuous and 
have continuous derivatives everywhere. In figure 7 we give an example of both flow 
geometries. 
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(a) 



(b) 



Figure 7: Representation of a ring of rolls with four vortex pairs (n = 4) and dimension- 
less annular width 5 = 0.5. Plot (a) shows radial profile of the dimensionless vertical 
velocity in a vortex for first flow (FW1, dashed line) and second flow (FW2, solid line), 
while plot (b) shows streamlines of horizontal fluid motion for both flows. Solid and 
dashed lines correspond to opposite circulations. 

Appendix B: Determination of S 

We consider an electromotive force £ generated by a flow structured in helical columns 
between two concentric cylinders that was coined a "ring of rolls". Assuming that u and 
B do not depend on z, we look for representations of £ in the general form 

roo 

S K (s)= / K KX (s,s')B x (s')s'ds', (B.l) 
Jo 

where k and A stand for s, if or z. Using a Taylor expansion of B x , we write the last 
equation as 

1 fYR 

S K (s) = a KX (s)B x (s)+(3 KXs (s)^-^(s)+--- (B.2) 

with 

roo 

a KX (s) = / K KX (s,s')s'ds', (B.3) 
Jo 

roo 

(3 kXs (s) = R K KX (s,s')(s'-s)s'ds'. (B.4) 
Jo 

The first term on the r.h.s of (B.2) represents the a effect, the second term represents 
the (3 effect, which will be omitted throughout the paper. The kernel K kX (s, s') depends 
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only on u. Under the first order smoothing approximation (FOSA), and using a definition 
of mean-fields by ip averaging, an analytical expression of K kX (s, s') was found in Avalos 
et al. (2007). The results are: 

2K ss (s,s) = — —^^{s,s)u v {s)u z {s) + -^{s,s)u z (s)u v3 {s)\, 



or/ / /\ f h n (s, s )~/\~//n.. h n (s, s ) „ . . „ , t . \ 
2K w (s,s) = — n u z {s)u s (s)^ u s (s)u z (s)\, 

1] \ S S J 



K sip = K^ s = K zs = K zz = 0. (B.5) 

The coefficients K sz and K vz are not zero, but the integrals / °° K sz (s, s') s' ds' and 
J °° K (pz (s, s') s' ds' can be shown to vanish. 
The Green's function h n are defined by 



= ^(7) ^ s'<s 



hn(s, s') 



1 (s 



2n \s 



for s < s' , 



(B.6) 



As in Avalos a/. (2007), we can further represent 



R 



m_l_ 



a KA if («A) = < 



(ss), (ipip) 



, q; k a = otherwise. 



(B.7) 

where a K \ is a dimensionless quantity independent of magnetic Reynold numbers R m ± = 
uo±R/v and R m \\ = uq\\R/t]. In figure 8, the s/R profile of the three non-zero dimen- 
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sionless a K \ coefficients are represented for both flows FW1 and FW2. 
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Figure 8: Radial profile of the dimensionless quantity a K \ for (a) FW1 and (b) FW2, both 
for n = 4 and 5 = 0.5. 



Appendix C: Numerical approach 

The correct handling of the non-local boundary conditions for the magnetic field is a 
notorious problem for the simulation of dynamos in non-spherical domains. Here, the 
kinematic eigenvalue problem in finite cylinders is solved by the integral equation ap- 
proach (Stefani et al. 2000, Xu et al. 2004a, Xu et al. 2004b, Xu et al. 2006). Basically, 
we use the following three integral equations: 



a o B(r')) x (r — r') , fiap r A(r') x (r — r') 



47r Jv |r - r'l 3 1 4ir Jv 



dV 



-^L mniC)x \hS dS '> (ai) 



|r-C 



l_ m = JL f (aoB(rQ).(C-O dr _jL/A(0.(C-0 



= 1 l B(r-)x(r Q 1 , B(0 

w 4vr Jy r - r' 3 47T 7s v y r - C v 7 
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where B is the magnetic field, A the vector potential, the electric potential, n the 
outward directed unit vector at the boundary S. The complex constant p contains as its 
real part the growth rate and as its the imaginary part the frequency of the eigenfield. The 
matrix a represents the a-effect defined by (5) or by (6). 

The reduction of the problem to cylindrical problems with azimuthal waves exp (i mip) 
was described in Xu et al. (2006). Finally, we end up with a generalised eigenvalue prob- 
lem for the critical dynamo number C c a (in the steady case), or for the complex constant 
p (in the unsteady case). The QR method is employed to solve this eigenvalue problem 
which gives also the eigenmodes of the magnetic field. 
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